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We have measured the intricate temperature dependence of the Co £2,3 x-ray absorption spectra 
(2p-3d excitations) of CoO. To allow for accurate total electron yield measurements, the material 
has been grown in thin film form on a metallic substrate in order to avoid charging problems usually 
encountered during electron spectroscopic studies on bulk CoO samples. The changes in spectra 
due to temperature are in good agreement with detailed ligand-field calculations indicating that 
these changes are mostly due to thermal population of closely lying excited states, originating from 
degenerate t2 g levels lifted by the spin-orbit coupling. Magnetic coupling in the ordered phase, 
modeled as a mean- field exchange field, mixes in excited states inducing a tetragonal charge density. 
The spin-orbit coupling induced splitting of the low energy states results in a non-trivial temperature 
dependence for the magnetic susceptibility. 

PACS numbers: 71.70.Ej, 75.10.Dg, 78.70.Dm 



It is an observation that many insulating 3d transition 
metal compounds do not reveal orbital degrees of free- 
dom in their electronic structure despite the incomplete 
filling of the 3d shell. This is related to the fact that the 
symmetry of the crystal often becomes so low that the 
degeneracy of orbitals is effectively lifted. This could be 
seen as a consequence of the Jahn- Teller theorem stat- 
ing that: "when the orbital state of an ion is degenerate 
for symmetry reasons, the ligands will experience forces 
distorting the nuclear framework until the ion assumes a 
configuration both of lower symmetry and of lower en- 
ergy, thereby resolving the degeneracy" i£ Nevertheless, 
for some compounds, with CoO and FeO as well known 
examples, a (nearly) cubic structure exists although these 
systems have a partially filled t2 g sub-shell. In such cases 
the orbital degeneracy can be lifted by the spin-orbit cou- 
pling (SOC). This in turn leads to quite interesting local 
physics, not only because of the formation of orbital mo- 
ments and large magnetocrystalline anisotropics, but also 
because of a strong temperature dependence in the local 
magnetic properties as we will address below. 

Taking the case of CoO, we will describe its ground 
state and low-energy local excitations using the theoreti- 
cal frame work laid out by Kanamori- and Goodenough 4 . 
The Co ions in this material are divalent and have 7 elec- 
trons in the 3d shell. Due to electron-electron interac- 
tions, a Hunds-rule S = 3/2 high-spin ground-state is re- 
alized, with an electronic configuration of approximately 
t\ g e 2 g ^ The single hole in the t2 g shell can be in three 
different orbitals, normally denoted as an xy, xz, or yz 
orbital. With the SOC being present, one should consider 
making complex linear combinations of these orbitals like 
^l/2(—xz — lyz) that carry an orbital momentum of 1 
Hb in the z direction. 6 Orbital momentum in the x or 
y direction can be created by cyclic permutation of the 



coordinates. One can therefore assign to each t2 g elec- 
tron a pseudo orbital momentum of / = 1 (note that it 
is still a d electron so I = 2.)2a4»£ For CoO this results in 
a total pseudo orbital momentum of L = 1. This orbital 
momentum couples with the spin S = 3/2 to a state with 
J = 1/2, 3/2, and 5/2. The doublet with J = 1/2 is the 
ground-state, with the quartet about 40 meV higher and 
the sextet about 120 meV higher than the ground- state. 
Furthermore, the sextet is split by the cubic crystal field 
into a quartet and doublet. In estimating the numbers 
we have used a SOC constant of £ = 66 meV from atomic 
Hartree-Fock calculations. 8,9 

Experimentally one could measure these excitations 
with the use of inelastic neutron scattering. The in- 
terpretation is somewhat complicated as one measures 
phonons, magnons and spin-orbit excitations all at the 
same time, but satisfying results have been obtained ] 10 ! 11 
Yet, it would be welcome to have an alternative spectro- 
scopic method with perhaps also advantages in terms of 
element specificity in case one is dealing with a multi- 
component materials, or higher material sensitivity for, 
e.g., thin films. In this regard it has been recognized that 
soft x-ray absorption spectroscopy (XAS) should have the 
potential to reveal the presence of those excited states: 
Tanaka and Jo ^ and De Groo^ have calculated that 
theoretically one should see a strong temperature depen- 
dence in the XAS lineshape. In such an XAS process, 
the soft-x-ray photon (770-800 eV) is used to promote 
a 2p Co core electron into the 3d valence shell. This 
transition is subject to strict dipole selection rules, e.g. 
A J = 0, ±1 in spherical symmetry, with the result that 
each particular initial state has its own set of reachable fi- 
nal states. It turned out that the differences between the 
sets of final states are appreciable when the differences in 
the initial states can be traced back to differences in the 
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FIG. 1: (color online) Theoretical isotropic Co-L2,3 XAS spec- 
tra for the different initial states split by spin-orbit coupling 
as present in CoO. 



quantum number J. In cubic symmetry as for CoO, J 
is not a good quantum number, but following analogous 
arguments each initial state with a different J will have 
a different spectrum. 

In Fig. 1 we have plotted the isotropic Co-L2,3 XAS 
spectra calculated for each initial state having a different 
J. Here we used the successful configuration interaction 
model that includes the full atomic multiplet theory and 
the hybridization with the O 2p ligands . 13 i 14 We have car- 
ried out the calculations for the Co 2+ ion in the Co06 
octahedral cluster using the XTLS8.3 program 14 with pa- 
rameter values typical for a Co 2+ system* 15 The spectra 
are dominated by the Co 2p core-hole spin-orbit coupling 
which splits them roughly in two parts, namely the L3 
(hv « 779 eV) and L2 (hv ~ 794 eV) white lines re- 
gions. The curves in Fig. 1 reproduce the main results 
of Tanaka and Jo, 1 — and De Groot We have not sepa- 
rately shown the two spectra for the doublet and quartet 
states originating from the J = 5/2 but combined them. 
The energy splitting between these two set of states is 
small enough compared to the energy difference between 
the ground-state and the center of the J = 5/2 states so 
that this splitting is of no real importance if the states 
are to be populated by temperature. These calculations 
show that the differences are so large that an experimen- 
tal resolution of order several 100 meV would be more 
than sufficient to discriminate the various initial states 
differing 100 meV or less in energy. Important is only 
the temperature and the statistics with which the spec- 



tra have to be measured. 

Nevertheless, it is rather surprising that to the best of 
our knowledge no experimental data have been reported 
in the literature so far to confirm the validity of the cal- 
culations, i.e. 14-17 years after the appearance of Tanaka 
and Joi 2 - and De Groot 1 ^ studies. Indeed, it turned out to 
be difficult to obtain reliable XAS data of CoO. Standard 
methods to measure the XAS signal failed so far: the to- 
tal electron yield (TEY) method cannot be applied since 
CoO is strongly charging at low temperatures, and the 
fluorescence yield (FY) method suffers from too strong 
self-absorption effects^ causing severe distortions in the 
spectral line shapes. We therefore have set out to pre- 
pare CoO in thin film form on metallic substrates, thin 
enough as to avoid charging problems upon measurement 
using the TEY method and yet thick enough to essen- 
tially retain the properties of bulk CoO. The thin film 
must also be polycrystalline, as to avoid complications 
caused by the occurrence of (magnetic) linear dichroism 
effects,Hl&l& 

The XAS measurements were performed at the Dragon 
beamline of the NSRRC in Taiwan. The spectra were 
recorded using the total electron yield method in a cham- 
ber with a base pressure of 3xl0 -10 mbar. The photon 
energy resolution at the Co £2,3 edges (hv « 770 — 800 
eV) was set at 0.3 eV. The actual polycrystalline CoO 
thin film was grown in-situ on a polycrystalline Ag sam- 
ple by means of molecular beam epitaxy (MBE), i.e. 
evaporating elemental and Co from alumina crucibles in 
a pure oxygen atmosphere of 10 _T to 10 -6 mbar. The 
base pressure of the MBE system is in the low 10 -10 mbar 
range. The thickness of the film was about 90 A as de- 
termined using a quartz-balance monitor, which in turn 
was calibrated using the time period of the oscillation in 
the intensity of reflection high energy electron diffraction 
(RHEED) pattern of CoO thin films epitaxially grown 
on Ag(001) j 19 i 2Q From that single crystal work we also 
know that a 90 A CoO film has a Neel temperature of 
about 290 K, i.e. essentially identical to the 291 K value 
determined for the bulk material. 22 The polycrystalline 
CoO thin film has been measured at normal incidence 
with a beam-spot size of about 1 mm 2 . Homogeneity of 
the thin film has been verified by measuring at several 
different positions on the thin film resulting in negligible 
differences of the spectra. 

Fig. 2 shows the Co-L2,3 XAS spectra of the poly- 
crystalline CoO thin film taken at various temperatures 
between 77 K and 400 K. The general line shape of the 
spectra shows the characteristic features very similar to 
that of bulk CoOi^, verifying the good quality of our 
CoO film. The low-temperature spectra have been mea- 
sured before and after heating to 400 K showing no de- 
tectable changes demonstrating the chemical stability of 
the sample. Important is that high quality spectra can 
be obtained for all temperatures, and this without the 
slightest indication for a reduction of the intensity for 
the lower temperature spectra, demonstrating that our 
thin film approach is a successful resolution against the 
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FIG. 2: (color online) Experimental Co-L2,3 XAS spectra of 
polycrystalline CoO thin film taken at various temperatures 
between 77 K and 400K. 
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FIG. 3: (color online) Theoretical Co-L2,3 XAS spectra of 
CoO at several temperatures. Temperature dependence is due 
to population of excited states according to Bolzman statis- 
tics. Non-cubic distortions and magnetic correlations have 
been neglected. 



charging problems often encountered when working with 
CoO. 

A closer look at the experimental spectra in Fig. 2 
reveals that there is a small but systematic change in 
the line shape with temperature. A blow-up of the spec- 
tra in the L2 edge range makes this even clearer. To 
understand this temperature dependence quantitatively, 
we simulated the spectra by summing up the theoreti- 
cal spectra starting from the different excited states (see 
Fig. 1), weighted by their thermal population according 
to the Boltzman statistics. These theoretical results are 
plotted in Fig. 3, and one can observe a very good agree- 
ment with the experiment. The predictions by Tanaka 
and Joi^ and De Grooti^ are thus confirmed. 

The calculations in Fig. 3 are done in cubic symme- 
try neglecting magnetic ordering. It is well known that 
CoO undergoes an anti- ferromagnetic ordering at 291 K. 
This phase transition is accompanied by a tetragonal dis- 
tortion and smaller distortions to even lower symmetry. 



Tanaka and Joi^ have shown that the non-cubic distor- 
tions (and the exchange field) can mix the low-energy 
states with different t^g orbital occupation, split by spin- 
orbit coupling, thereby changing the lineshape of the 
spectra. As the magnetic correlations and the tetrago- 
nal distortion are temperature dependent this could also 
result in a temperature dependence of the spectra. We 
will show below that the inclusion of these effects will 
even improve the already very good agreement between 
the temperature dependent theoretical and experimental 
Co £2,3 spectra as shown in Fig. 3. 

In the top panels of Fig. 4 we show the energy level 
diagram of the lowest 12 eigen-states of a CoOg 0- clus- 
ter as a function of exchange field (left) and tetragonal 
splitting between the d xy and d xz j yz orbital (right). The 
inset shows the ti g hole density. For cubic symmetry 
the ground-state is a Kramer's-doublet, which has a cu- 
bic charge density. For both an exchange field (in the 
z direction) and a tetragonal contraction of one Co-0 
bond, the d xy orbital becomes more occupied and the 
hole resides in a complex linear combination of the d xz 
and d yz orbital carrying an orbital momentum parallel 
to the spin of that hole. The bottom panels of Fig. 4 
show the changes of spectra (at K) as a function of 
these non cubic distortions. One can clearly see small 
changes in the spectra. The line-shape of the changes due 
to quantum mechanical mixing of the excited states into 
the ground-state are rather similar to the changes caused 
by thermal population. This leads to a compensating ef- 
fect: in going from low to high temperatures, one has a 
reduction of the exchange field and tetragonal splitting 
(they eventually disappear above the Neel temperature) , 
and thus changes in a direction opposite as those caused 
by the thermal population effect. The temperature de- 
pendence in the experimental spectrum is thus a result 
of the thermal population of excited states as shown in 
Fig. 3, slightly reduced in its magnitude by quantum 
mechanical mixing of the excited states into the ground- 
state at low temperatures as shown in Fig. 4. Indeed, a 
detailed comparison between Figs. 2 and 3 reveal that 
the theoretical spectra show too much change with tem- 
perature, and the inclusion of those quantum mechanical 
mixing effects from Fig. 4 will improve the agreement 
even more. 

The presence of the SOC has obviously important con- 
sequences for the ground state and the magnetic proper- 
ties of the Co ion. The SOC lifts the orbital degeneracy 
and the ground-state (without magnetic interactions) be- 
comes a Kramers doublet with an effective J = 1/2. One 
should note that this Kramers doublet is not Jahn- Teller 
active, that it has a cubic charge density and does not 
show linear dichroism. This Kramers doublet has a mag- 
netic moment, as revealed by the cluster calculations of 
L z -0.56a and S z -0.79ft (M z 2.15/ib, 9 4.29, 
and /ig ff = 3.7 2 /i#)- The induced Van Vleck moment is 
0.018/iB per Tesla. Exchange fields and crystal distor- 
tions mix in excited states, thereby altering the expec- 
tation value of the magnetic moment. At the same time 
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FIG. 4: (color online) Top panels energy level diagram of the 12 lowest states of a CoOq 0- cluster as a function of exchange field 
(left) and tetragonal crystal field distortion (right). Inset shows the ground-state t^ g hole density. Bottom panels: Theoretical 
Co-L2,3 XAS spectra of CoO as a function of exchange fields (left) and tetragonal crystal field distortion (right). 



these distortions induce a non-cubic local charge density 
and create linear dichroism. In order to calculate the low 
temperature ordered moments we need to make an es- 
timate for the exchange field and tetragonal distortion. 
CoO has a 1.2% tetragonal contraction at low T^i If 
we interpolate the values for the crystal-field and the ex- 
change field as found for CoO thin films grown on MnO 
and Ag 19 we find an exchange field of 12.6 meV and a 
tetragonal distortion of 25 meV between the d xy and the 
dxz/yz orbital. With these values we find for the ground- 
state that the total magnetic moments are S z = — 1.14ft 
and L z = -1.15ft (M z = 3.44/i B ). 

The SOC induced splitting of the low lying states has 
also important consequences for the temperature depen- 
dence of the magnetic susceptibility, in an analogous 
manner as it causes the temperature dependence of the 
XAS spectra. We have calculated the magnetic suscepti- 
bility x — hmjf^o dM/ dH (which is equal to Mj H for 
H small enough), for a CoOg 0- cluster. In Fig. 5a we 
show the the susceptibility calculated for the entire clus- 
ter including all Van Vleck moments (dotted line) and 
the susceptibility calculated for the lowest 12 eigenstates 
(straight line). For the second calculation we have thus 
neglected the Van Vleck (and Curie) susceptibility in- 
duced by the higher lying states, such as found around 1 
eV involving a t2 g -e g orbital excitation. This roughly cor- 
responds to experimentally subtracting a constant from 
the measured susceptibility such that the high temper- 
ature limit goes to zero. Panel b of Fig. 5 shows 1/x 
against the temperature. One can observe that the curve 



does not follow a linear behavior in a large temperature 
regime. There are two limiting cases of interest. For 
/c#T much smaller then the energy of the first excited 
state one should expect normal Curie- Weiss behaviour. 
This asymptote is indicated in Fig. 5b by a straight line 
on the left (blue). The other limit would be when ksT is 
much larger then the highest energy of the excited states 
split by SOC. In this limit one would again expect to re- 
trieve a 're-normalized' Curie- Weiss like behavior. This 
limit is indicated in Fig. 5b by a straight line on the 
right hand side (red). One can see that for temperatures 
starting roughly at 600 K 1/x shows linear behaviour and 
the high temperature limit is reasonably achieved. One 
might naively think that the high temperature region can 
be used to extract magnetic quantum numbers as it is of- 
ten done when a Curie- Weiss-like behavior is observed, 
but then one must also realize that those numbers are 
not representative for the ground state. 

To illustrate this more clearly, we will discuss both 
limits in more detail. We start with the low temperature 
limit. In Fig. 5c we plot /i^ ff , the (apparent) effective 
magnetic moment squared (/i^ ff is defined here as 3ks 
divided by the temperature derivative of l/x(T)). In 
the low temperature limit /i^ ff does not become a con- 
stant as one should expect for Curie- Weiss behaviour 
as we did not subtract the Van Vleck moments for the 
ground- state doublet. (The straight line is calculated 
for a subtraction of the Van Vleck moments for the low- 
est 12 states, which still allows the ground- state dou- 
blet to have a positive Van Vleck moment and the ex- 
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FIG. 5: (color online) Calculated magnetic susceptibility and 
its derivatives for a single Co 2+ ion in cubic ligand field sym- 
metry. Dotted lines without further correction, Straight lines 
Van Vleck (and Curie) moments induced by high energy (1 
eV and above) excited states neglected. The vertical lines 
mark the energies of the excited states. The horizontal lines 
mark the low (blue) and high (red) temperature limit of the 
Van Vleck subtracted susceptibility. Panel a: Susceptibility. 
Panel b: Inverse susceptibility. Panel c: The (apparent) effec- 
tive moment /^ ff . Panel d: The (apparent) Weiss temperature 
B as defined in the text. 



cited states a negative Van Vleck moment such that 
the sum of the lowest twelve states equals zero.) The 
linear behaviour of /i 2 ff for small T is easy to under- 
stand if one realizes that x(T) = /i 2 /(3/c#T) + Xvv, 
with xvv the constant Van Vleck susceptibility and 
fii the effective moment for the ground state. Then 
1/xiT) = (3k B T)/^ 2 - (3 2 k 2 B T 2 X vv)/^ + 0(T 3 ) and 
A4 = 3k B /d T (l/x) = Mi + 6xvvk B T + 0(T 2 ). 

The high temperature limit can be understood if 
one neglects Van Vleck contributions for all magnetic 
moments. One then finds that the average effective 



moment squared is just the average of the effective 
moments squared of each excited stated, weighted by 

their multiplicity: /i 2 N = Yli=\ v i^\ /J2iLi u i- Here Mat 
is the high temperature limit for the average effective 
moment squared of the first N non-degenerate sets of 
eigenstates, /a 2 the effective moment squared of the i-th 
non-degenerate set of eigenstates and Vi the degeneracy 
of the i-th set of eigenstates. For a cluster calculation 
of Co 2+ one finds a high temperature effective moment 
of /i 2 ff = AA 2 ji 2 B while the ground state value is 3.7 2 fi B . 
One does not only find a deviation from the ground-state 
effective moment but the high temperature limit also 
induces an apparent Weiss temperature. The Weiss tem- 
perature O is defined here as the intercept of the tangent 
to the l/x(T) curve with the abscissa. This can be calcu- 
lated by making a series expansion ofl/xinl/T. On = 

1 / k B (E*=i ^ A ;/E*=i "i ~ E*=i ^ A ^?/E*=i Wi), 
with Ai the energy of the i-th set of eigenstates. In panel 
d of Fig. 5 we plot the apparent Weiss temperature as a 
function of T. One can see a strong temperature depen- 
dence and the Weiss temperature slowly saturates at 242 
K. This apparent Weiss temperature has nothing to do 
with magnetic correlations since they were not included 
in this single ion calculation. Instead, the strongly 
varying (apparent) Weiss temperature merely reflects 
the fact that there are excited states with different 
magnetic quantum numbers thermally populated. 

Interestingly one finds that the high temperature limit 
for the (apparent) effective moment squared is reached 
at temperatures much lower then the energy of the high- 
est excited state. We have added vertical bars at the 
energies of the excited states in Fig. 5 and one can see 
that /i 2 ff becomes reasonable constant at a temperature 
that is about half the energy of the highest excited state. 
Not only does one find that the high temperature regime 
for /i 2 ff starts at rather low temperatures, one also finds 
that the low temperature region, where the ground-state 
Curie- Weiss and Van Vleck moments can be obtained, 
extends only up to 50K, roughly a tenth of the energy 
of the first excited state. This is particular disturbing as 
one naively would expect that excited states would only 
distort the susceptibility at temperatures comparable to 
the energy of the excited state. 

In the presence of magnetic correlations, such as in 
the real CoO material, those thermal population effects 
become very much masked or hidden. One then easily 
interprets the high temperature Curie- Weiss limit as a 
region in which one is free from magnetic correlations 
and yet still can capture the ground state local quantum 
numbers, and one also easily thinks that the measured 
Weiss temperature is only due to the magnetic correla- 
tions. The application of the Curie- Weiss law obviously 
leads to errors in case SOC is active. 

The effect of SOC excited states on the magnetic sus- 
ceptibility are well recognized in systems with an open 
4/ shells For 3d systems with an open t^g shell in nearly 
cubic symmetry these effects should also be important. 
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In that respect high-spin Fe 2+ , Co 3+ and Co 2+ are good 
candidates to show such effects. 

To conclude, we have been able to accurately mea- 
sure the temperature dependence of the Co £2,3 x-ray 
absorption spectra of CoO. The use of polycrystalline 
thin films of CoO grown on a Ag substrate allowed us to 
avoid charging problems usually encountered in electron 
spectroscopic studies on CoO. The changes in spectra 



as a function of temperature can be well explained in 
terms of the thermal population of closely lying excited 
states, originating from degenerate t<i g levels lifted by the 
spin-orbit coupling. The existence of these excited states 
also lead to a non-trivial temperature dependence for the 
magnetic susceptibility, e.g. one cannot apply the Curie- 
Weiss law to extract the relevant ground state quantum 
numbers from the high temperature region. 



1 H. A. Jahn and E. Teller, Proc. Roy. Soc. A161, 220 
(1937). 

2 H. A. Jahn, Proc. Roy. Soc. A164, 117 (1938). 

3 J. Kanamori, Prog. Theor. Phys. 17, 177 (1975); ibid. 17, 
197 (1975). 

4 J. B. Goodenough, Phys. Rev. 171, 466 (1968). 

5 S. Sugano, Y. Tanabe, and H. Kamimura, Multiplets of 
Transition- Metal Ions in Crystals (Academic, New York, 
1970). 

6 T. Jo, T. Shishidou, J. Phys. Soc. Jpn., 67, 2505 (1998). 

7 A. Abragam and B. Bleaney, Electronic paramagnetic res- 
onance of transition ions (Clarendon Press, Oxford, 1970). 

8 Given are the Hartree-Fock values for the spin-orbit cou- 
pling constant calculated with the program by R. C. 
Cowan. 9 £ is defined as the single particle coupling con- 
stant i.e. Hs.o. = J2i Ch ' s i- 

9 R. D. Cowan, The theory of atomic structure and spectra 
(University of Califormia press, Berkely, 1981). 

10 J. Sakurai, W. J. L. Buyers, R. A. Cowley, and G. Dolling, 
Phys. Rev. 167, 510 (1968). 

11 K. Tomiyasu and S. Itoh, J. Phys. Soc. Jpn. 75, 084708 
(2006). 

12 A. Tanaka and T. Jo, J. Phys. Soc. Japan 61, 2040 (1992). 

13 F. M. F. de Groot, J. Electron Spectrosc. Relat. Phenom. 
67, 529 (1994). 

14 A. Tanaka and T. Jo, J. Phys. Soc. Jpn. 63, 2788 (1994). 

15 Parameters for C0O6 cluster: A=6.5, £7^=6.5, £/cd=8.2, 
pda=-1.2, T pp =0.7, 10D<?=0.525, C=0.066; Slater integrals 
80% of Hartree Fock values; all values in eV. 



16 E. Pellegrin, N. Niicker, J. Fink, S. L. Molodtsov, A. 
Gutierrez, E. Navas, O. Strebel, Z. Hu, M. Domke, G. 
Kaindl, S. Uchida, Y. Nakamura, J. Markl, M. Klauda, G. 
Saemann-Ischenko, A. Krol, J. L. Peng, Z. Y. Li, and R. 
L. Greene, Phys. Rev. B 47, 3354 (1993). 

17 D. Alders, J. Vogel, C. Levelut, S. D. Peacor, T. Hibma, 
M. Sacchi, L. H. Tjeng, C. T. Chen, G. van der Laan, 

B. T. Thole, and G. A. Sawatzky, Europhys. Lett. 32, 259 
(1995). 

18 D. Alders, L. H. Tjeng, F. C. Voogt, T. Hibma, 
G. A. Sawatzky, J. Vogel, M. Sacchi, S. Iacobucci, and 

C. T. Chen, Phys. Rev. B 57, 11623 (1998). 

19 S. I. Csiszar, M. W. Haverkort, Z. Hu, A. Tanaka, H. H. 
Hsieh, H. -J. Lin, C. T. Chen, T. Hibma, and L. H. Tjeng, 
Phys. Rev. Lett. 95, 187205 (2005). 

20 S. I. Csiszar, thesis, U niversity of Groningen (2005); 
|http://irs.ub.rug.nl/ppn/287545644| 

21 W. L. Roth, Phys. Rev. 110, 772 (1958). 

22 J. Singer, Phys. Rev. 104, 929 (1956). 

23 Many examples of excited states due to spin orbit coupling 
effecting the magnetic susceptibility in 4/ compounds can 
be found. For example see: N. H. van Dijk, B. Fak, T. 
Charvolin, P. Lejay, and J. M. Mignot, Phys. Rev. B 61, 
8922 (2000). N. Metoki, K. Kaneko, T. D. Matsuda, A. 
Galatanu, T. Takeuchi, S. Hashimoto, T. Ueda, R. Settai, 
Y. Onuki, and N. Bernhoeft, J. Phys.: Condens. Matter 
16, 207 (2004). Y. Okada, Y. Inada, A. Galatanu, E. Ya- 
mamoto, R. Settai, and Y. O. Nuki, JPSJ 72, 2692 (2003). 



